Anderson Impurity in Helical Metal 
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We use a trial wavefunction to study the spin-1/2 Kondo effect of a helical metal on the surface of 
a 3-dimensional topological insulator. While the impurity spin is quenched by conduction electrons, 
the spin-spin correlation of the conduction electron and impurity is strongly anisotropic in both 
spin and spatial spaces. As a result of strong spin-orbit coupling, the out-of-plane component of 
the impurity spin is found to be fully screened by the orbital angular momentum of the conduction 
electrons. 



On the surface of a 3-dimensional (3D) topological in- 
sulator, massless helical Dirac fermions emerge, [lj Spe- 
cific materials of 3D topological insulators have recently 
been studied and observed. |2H5| Interesting properties of 
these helical Dirac fermions have become a focus of re- 
cent studies. In a helical metal, spins are coupled to 
momenta, so that magnetic properties are expected to 
be highly non-trivial. Theoretical study on this aspect, 
however, has so far been limited to the effect of a classical 
magnetic impurity. [10( 
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The low temperature property of a quantum spin-1/2 
magnetic impurity or an Anderson impurity in a con- 
ventional metal|ll| has been an interesting and impor- 
tant problem for decades in condensed matter physics. 
The problem has been studied by using the renormal- 
ization grou p Il2 | . the Bethe ansatz[l3l [Tij ]. the 1/iV- 
expansionjla. Il6j . and the conformal field theory |l7|. 
The phenomenon is a well known Kondo effect, in which 
the impurity spin is completely screened by the conduc- 
tion electrons. Because of the coupling between spins 
and momenta, the Kondo effect in a helical metal will be 
interesting to be examined. In this paper, we use a varia- 
tional method to address this problem. We show that the 
correlation of the impurity spin and the conduction elec- 
tron spin density is strongly anisotropic in both spatial 
and spin spaces, in contrast to the isotropic screening in a 
conventional metal. While the impurity spin is quenched 
at low temperatures, similar to the usual Kondo prob- 
lem, the correlation between the impurity spin and the 
spin of the conduction electrons in the helical metal is 
significantly reduced. The reduction of the spin-spin cor- 
relation in the direction perpendicular to the surface can 
be shown explicitly to be compensated by the screening 
of the orbital angular momentum of the conduction elec- 
trons. The possible experimental consequences will be 
briefly discussed. 



We consider a spin-1/2 magnetic impurity in a helical 



Figure 1: (color online), (a) Schematic energy- momentum 
dispersion in a helical metal. Electron spin (green arrow) is 
perpendicular to momentum, (b) Energy diagram of a heli- 
cal metal with a singly occupied Anderson impurity state at 
energy in the absence of their coupling H m i X in Eq. fTJ) . fi 
is the chemical potential. 



metal in a 2D x-y plane, described by the Hamiltonian, 



H — H c + H mix + Hd 

H c = cJJfivF {cr X k) • z- n]cu 



(1) 
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H d = (ed - n){ridf + Udx) + Un d ^n d \,- 

In the above equations, Ck and d are annihilation oper- 
ators of the conduction electron with momentum k and 
of the (i-electron at impurity site in the spinor represen- 
tation, respectively, cr are the Pauli matrices, \i is the 
chemical potential, and rido = d\d a . H c describes a heli- 
cal metal, and Hd is a local impurity Hamiltonian with 
the impurity energy level and U the Coulomb repulsion 
energy of the two d-electrons on the same impurity site. 
H m i x is a hybridization term between the helical metal 
and the impurity state. 

To start with, let us discuss the simple case H m i X = 
first. In this limit, the helical metal and the impurity 
state are decoupled. The single electron eigen-energy and 
its corresponding eigenstate of the helical metal H c are 
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given by 



e k± = ±hv F \k\ - fx 



1 



7k± = -j=(e^ k c kt ±ie~ 



(2) 



with # k the angle of the momentum k with respect to the 
x-axis that tan # k = k y /k x and "±" refer to upper and 
lower bands, respectively. For each single particle state, 
its spin lies in the plane and is perpendicular to the direc- 
tion of its momentum. Note that such a spin-momentum 
relation has recently been reported in a spin- and angle- 
resolved photoemission spectroscopy (spin-ARPES) ex- 
periment [l8|]. The ground state of H c is then given by 



n tLi°>. 
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where the product runs over all the states within the 
Fermi sea f2. As for the impurity part, we shall consider 
an interesting case, where td < A* < e d + U, so that the 
impurity site is singly occupied and has a local moment. 
The total energy of the system Hq = H c + Hd is then 
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where s = ± is the band index and hereafter the sum of 
{ks} is always over the Fermi sea. In Fig. Q] we illustrate 
the electron dispersion and the spin of the Dirac fermion 
and the ground state of a helical metal, in the absence of 
the hybridization with the impurity state. 

To study the ground state of H in the presence of H mix , 
we will use a trial wavefunction approach. Such a trial 
wavefunction method was used to study the ground state 
of the Anderson impurity problem in the conventional 
metal [l6l.lT9j or in an antiferromagnet [2(| . For simplicity, 
we consider here the large U limit, and exclude the double 
d-electron occupation at the impurity site. We expect the 
qualitative physics should apply to the case of finite but 
large U. The trial wavefunction for the ground state is, 



I*) = (ao+ ^a ks 4 s 7k s )|*o> 

{ks} 



(3) 



where d k ± — -^{ez 9k d^ ± ie ^ Sk di), and an and a ks 
are variational parameters . They are to be determined 
by optimizing the ground state energy. Note that by 
construction, the impurity state is either singly occupied 
or non-occupied. 

The energy of H in the variational state |\&) is given 

by 

12iks}( E o ~ £ ks)aL + 2V k a a ks + e ks ag 
E = — - — - s — — s (4) 
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The variational method leads to the equations below, 

e ks )a = ^kflks (5) 

{ks} {ks} 

(E - E + e ks )a ks = V k a (6) 



which may determine the binding energy A& = Eq — E 
due to the hybridization, 



E 

{ks} 
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If Ah > 0, the hybridized state is stable against the de- 
coupled state. To proceed further, we introduce an en- 
ergy cutoff A for the helical metallic state in our analysis, 
and consider Vk = VQ(A — hvp\k\), where <d(x) is a step 
function, which is 1 for x > and for x < 0. The low 
energy physics is expected to be insensitive to the cutoff 
A. A natural energy cutoff for the helical metallic state 
in the topological insulator is the half of the bulk gap. [|[ 
Eq. (JT]) enables us to determine the ground state energy E 
or the binding energy A;, and the ground state wavefunc- 
tion. The results are that the hybridization always leads 
to a binding A^ > 0, and a magnetic screening of the con- 
duction electrons to the impurity spin at any =^ 0. At 
the Dirac point /j, = 0, the binding and screening only oc- 
cur at the dimensionless hybridization strength a — jjris 
above a critical value. We will come back to this point 
later. 

We shall first examine the magnetic properties of the 
system, which are most interesting. The central quantity 
we will examine is the correlation function of the impurity 
spin Sd = ^cftard at the impurity site (set to be at the 
origin r = 0) and the conduction electron spin density 
S c = ic^(r)crc(r) in the ground state, namely 



J uv {v) = (S^(r)SS(O)), 



(8) 



where (Q) is the ground state average of Q, and u, v = 
x,y,z are the spin indices. If u and v are both on the 
x — y plane, we then have, by the rotational symmetry, 
J uv (r) = J„'„'(r'), if u' = R z {p)u, v' = R z ((3)v, and 
r' = R z (f3)r, with R z ((3) a rotational operator of angle (3 
along the z-axis. J^ yV (r) can be calculated by using the 
trial wave- function in Eq. (3). Note that the variational 
parameter a k ± is independent of # k since it can be ex- 
pressed as Vao/(±hvF\k\ — [i — A&) according to the Eq. 
©. In terms of a k ±, the diagonal components are found 
to be 

Jz,(t) = -I|£(r)| 2 + i|.4(r)| 2 
■Mr) = -l\B(r)\ 2 -±Re(A 2 (r)) 
J yy {v) = - l -\B{v)f+ l -Re{A\v)) 



where A(r) = £{ks} se^+^Vs and B(r) = 
D{ks} e lk r a ks . J 2Z (r) = J zz (r) is rotational symmet- 
ric around the impurity site, while ./^(r) and J vv (y) are 
highly anisotropic in space. In Fig. [2] we show J zz (r) 
and J yy (x,0). Note that J xx (x,0) — J zz (x,0) and their 
signs oscillate in space. A negative (positive) value means 
anti-parallel (parallel) correlation between the impurity 
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spin and the conduction electron spin density. Due to 
the absence of the spin SU (2) symmetry, the off-diagonal 
components are not zero in general and they are 



Jxy{?) — Jyx{ r ) 



J xz (r) = -J«(r) = -Im(B(r)A(r)), 



J yz {r) = -J zy (r) 



-Re{B{v)A{v)). 
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Figure 2: (color online). Spatial Spin-spin correlation func- 
tions J yy (x,0) and J zz (x,0) (J xx (x,0) = J zz (x,0)) in Eq. (8) 
as functions of x, for a set of parameters Ab = 0.005A and 
fi = — 0.01A. The length unit is Iivf/A. Different parameters 
show qualitatively similar feature as long as A;, > 0. 

The spatial distributions of the spin correlations in 2D 
are plotted in Fig. [3j As we can see, the magnitude of 
J yy (r) is larger at the space point r near the x-axis than 
near the y-axis, and the magnitude of J xx {r) is larger at 
the space point r near the y-axis than near the x-axis. 

The anisotropic spin correlation function implies that 
the Kondo screening to a magnetic impurity in a helical 
metal has more complex texture. In the strong spin-orbit 
coupled system, the screening to the impuriy spin may 
be contributed from both the spin and orbital angular 
momentum of the conduction electrons. In what follows, 
we will show that in comparison with the Kondo effect 
in a conventional metal, the spin-spin screening is largly 
reduced, and the z-component of the impurity spin is 
fully screened by the orbital angular momentum. Let us 
first examine the correlation between the total spin of the 
conduction electron, S c = JS c (r)d 2 r and the impurity 
spin. We define 



(S u c S u d )/ 



{ks} 



with u = x,y,z, where the denominator is the occupa- 
tion number of the d-electron. / = I u is a measure 



Figure 3: (color online). Spatial spin-spin correlation 
functions plotted in x-y plane, (a) J zz (x, y); (b)J xx (x,y); 
{c)J yy (x,y); (d)J xy (x,y) or J yx (x,y); (e)J xz (x,y) or 
— Jzx{x, y)); {f)J V z{x,y) or — J zy (x,y). The parameters are 
the same as those in Fig. [2] 



of spin-spin screening strength. J = if the conduction 
electron and impurity state are decoupled. I = —3/4 for 
a spin-f/2 magnetic impurity in a conventional metal, 
corresponding to a spin singlet of two spin-1/2. In the 
present helical metal, when the Fermi level is at or below 



the Dirac point, we find I x = I y 



■i and I z = 0, there- 



fore I = —j, which is one third of the value of —3/4 in the 
usual Kondo problem. Note that although J zz is non-zero 
locally, its overall contribution to the spin-spin screening 
is zero. The spin-spin screening comes from the in-plane 
spin components. Since / is independent of any parame- 
ters including the energy cutoff (the only requirement is 
fj, < 0), we expect that this is a general property of the 
helical metal. When the Fermi level is above the Dirac 
point, there is a deviation from the above result. For in- 
stant, with the parameters /i = 0.01 A and = 0.005A 
the deviation is about sixty percents of — f/4. At \i > 0, 
the upper Dirac cone is partially filled. The electron state 
in the upper Dirac cone has spin orientation opposite to 
the corresponding state in the lower Dirac cone, which 
provides a stronger screening than induced by the lower 
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band only. We remark that the difference between fi < 
and \i > is due to the particle-hole asymmetry in our 
model parameter with U large. If U + 2(e d — /x) = 0, we 
have particle-hole symmetry, and we expect correponding 
symmetries for fi < and \x > 0. 

The conduction electron spin is not a good quantum 
number in a helical metal. To better understand the 
Kondo screening in the helical metal, we shall consider 
contribution from the orbital angular momentum. Be- 
cause of the two dimensionality, only z-component of the 
orbital angular momentum can be considered here. The 
system is invariant with respect to a simultaneous rota- 
tion of both spin and space along the z-axis at the origin 
of the impurity site. Therefore, the z-component of the 
total angular momentum J z = L z + S* + S d commutes 
with the Hamiltonian H and is a good quantum number, 
where L z is the total orbital angular moment of the con- 
duction electrons. Since the ground state preserves the 
time reversal symmetry, we have J z \^>) = 0. Because of 
this and (S z S d ) = discussed above, we have 

(L z S z d ) _{J z S z d -S z S z d -{S z d f) _ 1 



S{ks} a l 



(9) 



Therefore, although spin-spin screening in the z-direction 
in the present Kondo problem vanishes, the orbital an- 
gular momentum L z replaces the role of the conduction 
electron spin S z to screen the impurity spin. 

The spatial anisotropic correlations may be detected 
in spin resolved scanning tunneling spectroscope (STM) 
experiments. The magnetic susceptibility is finite in this 
variational theory. In the light of the unconventional 
spin correlations, exactly how the susceptibility behaves 
has to be answered by more elaborate treatments such 
as quantum Monte Carlo and numerical renormalization 
group. This issue is left for further investigations. 

We now discuss the binding energy of the system. 
From Eq. (JT)), we obtain 



(d - ^ 

a(ji + A b ) In 



A 6 + a(A- \n\) 

4(A + /i + A 6 )A 6 



(ji + |m| + 2A fc ) 2 
In the limit of small a, and a A < /i — e d , we have 



(10) 



At ~ A exp 

2 

A fc « £jr- exp 



fj, — eg — aA 



,/i < 

,/i > 0. 



(11) 



If the Fermi level is at the Dirac point, /i = 0, A& has a 
positive value solution only if a > a c — |ed|/A. In other 
words, for a given hybridization strength a, there is a 
critical value of the impurity d-level, below which there 
is no magnetic screening. This is due to the vanishing 
density of states at the Dirac point. This result is simi- 
lar to the impurity problem in graphene. pil. [22j| . The 
binding energy as functions of the chemical potential fi 
and the hybridization a are plotted in Fig. 21 
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Figure 4: (color online). Binding energy A;, as functions of 
the chemical potential p and the hybridization strength a. A 
is an energy cutoff. 



Ab is asymmetric about the point fi — 0. The solid 
curves in Fig. 0] are for a = a c and \x — respectively. 
At a > a c , the hybridized state is always stable. At 
a < a c , the binding energy is strongly reduced around 
the Dirac point and vanishes at fj, = 0. As a — > a c + + , 
Ai, ss \a — a c \A. This behavior is consistent with the 
result by using large-degeneracy method. [23| 

In summary, we have examined a spin- 1/2 Ander- 
son impurity in a 2D helical metal. The momentum- 
dependent spin orientation in the helical metal shows 
interesting magnetic properties. We have used a trial 
wavefunction method to study the system at the large 
Coulomb U limit. While the impurity spin is quenched by 
the conduction electrons, we find strong spin- and spatial- 
anisotropies in the correlations between the impurity spin 
and the conduction electron spin density. Because of the 
strong spin-orbit coupling, the orbital angular momen- 
tum also contributes significantly to the screening of the 
impurity spin. In particular, the out-of-plane component 
of the impurity spin is found to be fully screened by the 
orbital angular momentum of the conduction electrons. 
After our paper was submitted, there was a report on the 
experimental realization of the magnetic doping on the 
surface of the 3D topological insulator Bi2Se3 [2J|, and 
a e-print by Zitko, who studied the essentially the same 
model as ours and mapped it onto a conventional Ander- 
son pseudogap impurity model. [25| He showed that the 
impurity spin is fully screened. Our trial wave-function 
method reported here is consistent with his result. 

Finally we comment on the similarity and difference 
of the helical metal with the graphene. In terms of the 
eigen-energy problem, a helical metal is often considered 
to be a quarter of a graphene, where the pseudospin 
plays the similar role with the spin in helical metal. The 
graphene has a two-fold spin degeneracy and two-fold de- 
generacy associated with the two sublattices. Our results 
for the binding energy are similar to the graphene. How- 
ever, because of the difference between the spin in the 
helical metal and the pseudospin in graphene, the mag- 



5 



netic properties of the two are different. If we consider 
only one of the two Dirac cones is coupled to the spin- 
1/2 impurity, the graphene behaves like a conventional 
metal. 
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